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ABSTRACT 

Context. With the discovery during the past decade of a large number of extrasolar planets orbiting their parent stars at distances 
lower than 0. 1 astronomical unit (and the launch and the preparation of dedicated space missions such as CoRoT -Borde et al. , 2003- 
and KEPLER -Borucki et al, 2007-), with the position of inner natural satellites around giant planets in our Solar System and with 
the existence of very close but separated binary stars, tidal interaction has to be carefully studied. 

Aims. This interaction is usually studied with a ponctual approximation for the tidal perturber. The purpose of this paper is to examine 
the step beyond this traditional approach by considering the tidal perturber as an extended body. To achieve this, the gravitational 
interaction between two extended bodies and more precisely the interaction between mass multipole moments of their gravitational 
fields and the associated tidal phenomena are studied. 

Methods. Use of Cartesian Symmetric Trace Free tensors, of their relation with spherical harmonics and of the Kaula's transform 
enables us to analytically derive the tidal and mutual interaction potentials as well as the associated disturbing functions in extended 
bodies systems. 

Results. The tidal and mutual interaction potentials of two extended bodies are derived. In addition, the external gravitational potential 
of such tidally disturbed extended body is obtained, using the Love number theory, as well as the associated disturbing function. 
Finally, the dynamical evolution equations for such systems are given in their more general form without any linearization. We also 
give, under simplified assumption, a comparison of this formalism to the ponctual case. The non-ponctual terms have to be taken 
into account for strongly deformed perturbers (72 ^ 10^) in very close systems (ob/^b ^ 5). 

Conclusions. We show how to derive explicitely the dynamical equations for the gravitational and tidal interactions between extended 
bodies and associated dynamics. Conditions for application of such formalism are given. 

Key words. Methods: analytical - Celestial mechanics - Planetary systems - Stars: close binaries 



1. Introduction 

In celestial mechanics, one of the main approximation done in the modelling of tidal effects (star-star, star-innermost planet or 
planet-natural satellites interactions) is to consider the tidal perturber as a point mass body. However a large number of extrasolar 
Jupiter-like planets orbiting their parent stars at a distance lower than 0. 1 AU have been discovered during the past decade (Mayor 
et al, 2005). Moreover, in Solar System, Phobos around Mars and the inner natural satellites of Jupiter, Saturn, Uranus and 
Neptune are very close to their parent planets. In such cases, the ratio of the perturber mean radius to the distance between the 
center of mass of the bodies can be not any more negligeable compared to 1 (cf. Fig. [T]l. Furthermore, it can be also the case 
for very close but separated binary stars. In that situation, neglecting the extended character of the perturber have to be relaxed, 
so the tidal interaction between two extended bodies must be solved in a self-consistent way with taking into account the full 
gravitational potential of the extended perturber, generally expressed with some mass multipole moments, and then to consider 
their interaction with the tidally perturbed body. In the litterature, not so many studies have been done (Borderies, 1978-1980; 
Ilk, 1983; Borderies & Yoder, 1990; Hartmann, Soffel & Kioustelidis, 1994; Maciejewski, 1995). The purpose of this work is 
then to provide a theoretical procedure to obtain this tidal gravitational interaction as well as its associated tidal dynamical evolution. 

Several years ago, Hartmann, Soffel & Kioustelidis (1994) introduced in Celestial Mechanics an interesting tool, based on 
Cartesian Symmetric Trace Free (STF) tensors, to treat straighforwardly the couplings between the gravitational fields of extended 
bodies. These tensors are fully equivalent to usual spherical harmonics but in addition a set of STF tensors represents an iiTeductible 
basis of the rotation group S03 (Courant & Hilbert 1953, Gelfand et al. 1963). It means that using algebraic properties of STF 
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tensors with the index notation of Blanchet & Damour (1986), these objects become a powerful tool to determine the coupling 
between spherical harmonics in an elegant and compact way. However, as these tensors are not widely used in Celestial mechanics, 
we first recall their definition and fundamental properties and we stress their relation with usual spherical harmonics. Then, we treat 
the multipole expansion of gravitational-type fields for which each type is related to a given extended body. First, the well-known 
external field of such body is derived using STF tensors; classical identities are provided. Next, the mutual gravitational interaction 
between two extended bodies and the associated tidal interaction are derived. We show how the use of STF tensors leads to an 
analytical and compact treatment of the coupling of their gravitational fields. We deduce the general expressions of tidal and mutual 
interaction potentials expanded in spherical harmonics. Using the classical Kaula's transform (Kaula, 1962), we express them as a 
function of the Keplerian orbital elements of the body considered as the tidal perturber. These results are used to derive the external 
gravitational potential of such tidally perturbed extended body. Introducing a third body, its mutual interaction potential with the 
previous tidally perturbed extended body is defined that allows us to derive the disturbing function using the results obtained with 
STF tensors and the Kaula's transform. At this stage, the different type of mutual gravitational interaction are defined. The dynamical 
equations ruling the evolution of this system are obtained. Finally, we use a reduced form of that equations to quahtatively quantify 
the influence of non-ponctual terms of the disturbing function in comparison with the ponctual case. 




Fig. 1. System of two extended bodies. Ra and are the respective mean radius of A and B while Dab is the distance between 
their respective center of mass, sa and sb are the conical angles with which each perturber is seen from the center of mass of the 

perturbed body. When ^ 1 with / =A or B the i* body could be considered as a ponctual mass perturber. 



2. STF-Multipole expansion of gravitational potentials 

2.1. Definitions and notations 

We focus first on the STF tensors. Let us define a Cartesian 1-tensor as a set of numbers T,,,, with / different indices ii to //, each 
taking integer value running between 1 and 3. A compact multi-index notation, first introduced by Blanchet & Damour (1986), is 
generally used. An uppercase latin letter denotes a multi-index while the corresponding lowercase denotes its number of indices : 

L = /i/2...//, = r,-, . (1) 

The Einstein summation convention is assumed in the following, so if some index appears twice, a summation over that index is 
implied 

AkBk - Ai^i^_,ji^Bi^i^_,ji, - ^ Ai^i^_,,ji,Bi-^i-,,ji_ . (2) 

hJi,--Jk 

Given a Cartesian tensor we denote its symmetric part with parenthesis 

■^(L) = ^(h...ii) = ^ ^ '^V(i)iV(2)...io-<n ' (3) 

O" 

where cr runs over all the II permutations of { 1, 2, 3, /}. 
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Next, the Symmetric Trace Free part of a tensor is denoted indifferently by - yi<L> = -^oiij i,>- Following Thorne 
(1980), the STF part of Ji^ reads 

[|] 

— ^ ' ^(iik-'-^ht-lht^ ilk+i-'l)a]ai..atat (4) 



where 6 is the classical Kronecker delta function. 



Sl - ^(L) and a 



Kl 



II (-If (21 - 2k - ly.i 

(2l-l)\\ (2k)\\ (l-2k)\ 



j denoting the integer part of | while 

l\^l(l-l)(l-2)---2xl and /!! = / - 2) (/ - 4) ■ ■ ■ x 



1 if / is odd 

2 if Hs even 



(5) 



2.2. STF-basis 

Let e, (/ running between 1 and 3) be a Cartesian basis vectors set (with the rule e, ■ - 6ij). The basis of the (21 + 1 )-dimensional 
vector space of STF rank /-tensors is made of the STF parts of the /-fold tensorial products (Thorne, 1980) 



m 












(g)E» 


(6) 


- n=l 




_p=m+ 1 





where 



with f = -1. For m > 0, let us define the following algebraic object 



n^-: UK 



p=m+l 



Then, the STF canonical basis is proportional to £J2,> '-^'^ chosen as (Thorne, 1980): 



{//,m _ ilniplm 



where 



a'" = (-1)'"(2/- 1)!! 



21+1 



4n(l - my.(l + my. 



The constant A'"' is choosen to get a normalization such that 



ifi,m/ifi,,n'Y ^ (21 +m 



Anil 



(7) 



(8) 



(9) 



(10) 



(11) 



z* corresponds to the complex conjugate of z, z being a complex number or function. Finally, taking into account Eqs. Q, Q and 
( [TO) i, we obtain 

[¥] 



where 



.'- = .4=^.^?^ and M. = (-l)" 



21+ 1 (I -my 



2'j[(l-jy(l-m-2jy 



An (l + my' 



(13) 



Let us consider an arbitrary vector x = (x', x^, x^). We can now introduce the Euclidean norm r of x, the corresponding unit vector 
n and the Cartesian tensors x^ and n'^ constructed on x as follow 



r , x^ = x'' x'^...x" , = n'' n'^...n" 



Using the harmonic property = (for r > 0) and Eq. |l4| , one gets 

dJ-] = dLl-]^(-i)'(2i-iy 



„/+i ■ 



(14) 



(15) 
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where di and hi are the STF part of 5/. and n^, respectively. Then, considering Eq. |l2| i and taking into account ('"^) = 0> the 
computation of Kronecker deltas functions combined with Eq. ( [T5| ) leads to 



•^TM-j = N„„a""^d, + idyrd[ 

= (-iy(2/-l)!!^ 



By applying twice Eq. ( 16 1, we find 

'1 



Njkaj''^\d, + idyfdi 



This last equation leads to a composition law of a product of basis functions J'/". After some algebra, we obtain 



l.m 



■ LJ 



where 



-v''f = y^'' 

'jk ' Im 



21+1 



2j +1 [(/ + - (m + k)] ! [(/ + j) + (m + k)] ! 



; (I + m)\(l - m)\ ij + k)\{j - k)\ An [2 (I + j) + 1] 

Taking into account that the left hand side of the Eq. ( fTS) can be written using Eq. ( [T5| as follow 

'V 
rl 



= (2/ + 2; -l)!!i^;-'"i^f-^, 



we finally get 



v'.- , = (2^-l)'K2;-l)!! , 
■^L'^j^u (2/ + 2j-l)!! ^-z*^" 



(16) 



(17) 



(18) 



(19) 



(20) 



(21) 



2.3. Relation between STF-basis and usual spherical harmonics 

In this section the equivalence between STF basis and spherical harmonics Yi^„ is given. Following Abramowitz and Stegun (1970), 
Yi^m reads 



Yi miO, (fi) - 7V/,„P"'(cos ff) exp [im(p] for m > 

[¥] 

= M,,, (exp [iV] sine)"' ^ fl'"'-'(cos e')'-"'-^^ , 

where the P™ are the classical associated Legendre Polynomials. We also recall the symmetry property of the Vi m, namely: 

¥,,-,„ (e, if) - {-IT y,;, (0, if) with m > 0. 
Since the components of the unit vector n in complex form as follow 

rix + iriy - exp {i(f \ sin , n- - cos 6 , 

we can identify 



which can be inverted by using Eq. ( [TT| l as 



Anil 



+/ 



(2/+ 1)! 



m=-I 



To give an example of how it works, let us examine the cases / = and I - I. The spherical harmonics are 

1 

Yo,o 



Fin = 



Fi 1 = - ■ 



sin exp [iip] . 



(22) 

(23) 

(24) 
(25) 

(26) 

(27) 
(28) 
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Let us compare with the STF-basis functions. For I - 0, the normaUzation rule given in Eq. (Ill gives the single number 1 / V47r. 
For / = 1 we have 



1,0 \ 



, 2 



o^ 



1 J 



while 



J/ 

V 3 / 



1 ^ 
/ 

oj 



and 



(if 



1,-1 \ 



1,-1 



1 ^ 
— / 
j 



(29) 



it is verified that Y^m = S/j"\x'/r). 



2.4. Multipole expansion of the external gravitational field of an extended body 




<Pab 



Fig. 2. Spherical coordinates system associated to the equatorial reference frame 'Re : {(9a,Xe, Ye,Ze) of an extended body A; 
we have r = (r, 6, tp) and r^s = (fAB, Sab, Vab) where r^s, Oab and ^ab are the coordinates of the center of mass of the potential 
extended perturber B. /?a is the equatorial radius of A. 



Let us consider some matter distribution, corresponding to a body A, in the inertial coordinates (f, x'). The Newtonian gravita- 
tional potential of this body, V'^ (f , x), is obtained by solving the Poisson equation 

(f, x) = -47tGpa it, x) with lim V^{t, jc) = , (30) 

\x\-Kxy 

Pa (t,x) being its density. This leads an expression of |jt:| > Ra, where Ra is the equatorial radius of A, to 



cfP^d^x'. (31) 
Ja \x-x'\ 



Then, the external gravitational field of the body A for |r| > Ra can be represented by the series; 



Ia=0 

where all mass multipole moments are defined by: 



V\t, X) = G ^ ^-^MlA. - h (32) 



Ml. = n[^, t ^'^^-^^Lr (33) 

"ia=-'a 



the usual gravitational moments in the physical space being given by: 



An r 

Mi.,m. = r'^ Yi,„, (0, V) dMA; (34) 
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Ma is the mass of A and (JMa = PA'dr sin OdOdtp. Inserting Eq. ( 33 i into Eq. ( 32 1 and using Eq. ( 16 1, the final expression for V 
for |r| > Ra is thus obtained: 

Ia=0 mA=-lA 

One should note the symmetry property of Mi^^m^: 



yiA,mA (0,<P) 



(35) 



M,„_„, = (-ir-M,;,„^. 

Moreover, Mi^^,„^ could be represented in its polar form 

^/a,«u = \MiA,mJ exp [i6Mi^j„J 
where the following identities are obtained from Eq. (|36]l 



\M,A,-n,A\ = \MlA,mJ 

Arg (M/^ _„,^) = otatt - Arg (M/^,„,^) = matt - SMi^^,„^ 



(36) 



(37) 



(38) 



Using the classical symmetry property concerning spherical harmonics given in Eq. (23 i, (t, r) could also be expressed with the 
associated Legendre polynomials: 



" ^ p;'Mcos6») 
V (t,r) 2_j ^ [C/A,mA cos (mA<p) + 5/^,,,,^ sin (wa^)] : 

/a-0 mA-0 

where the usual coefficients C/^,m^ and S i^^m^^ are given by: 

C/A,mA = - <5„u,o) Re (M,„,„J 

= -2N';'^' (1 - 5„,„o) Im (M,,,„J ■ 

The expression of Mi^ ,„^ and 6Mi^ „,^ are then deduced for toa > 0: 



\MlA 



c. 



'a.wia 



(2 - 5,„A,o) 



> /A,mA 



(1 - ^mA.of 



6M, 



I A, Ma 



-Arctan 



(1 - <5mA,o) (2 - ^niA.o) 



/a./iia 



In the general case, the gravitational moments are expanded as: 



lA.WlA /a.jHa /a.BIa 



(39) 

(40) 

(41) 
(42) 

(43) 



Mf^^^ and M^*^^^ are respectively those in the case where A is isolated (without any perturber) and those induced by the tidal 
perturber(s). 

One can identify some special values of Mi^ „j^ relevant for the gravitational field of a body A. The trivial one is its mass. Ma 

Mo ft = ViTrMA. (44) 

Furthermore, we know that the external field of an axisymmetric body A can be expressed as a function of the usual multipole 
moment 7/jM(see e. g. Roxburgh 2001) 



V^(f,r) = 



GMa 



'a>0 



using Eq. ([35]), we identify in a straigthforward way: 



V'^{t,r)^G 



Ma 



y ^'a.o 



Ia=0 



where ., „ = Mf* + M^ 

■'/a.'a.U /a,0 'a 



K ■ 



(45) 



(46) 



We can now focus on the second type of gravitational interaction, namely the tides between two extended bodies. 



' They are driven by two types of deformation. The first one is those induced by internal dynamical processes such that rotation (through the 
centrifugal acceleration) and magnetic field (through the volumetric Lorentz force). The second one is the axisyrmnetric permanent tidal oval shape 
due to a companion in close binary or multiple systems. 
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2.5. Determination of tiie tidal potential 

Let us now introduce an accelerated reference frame, i.e. (t, X^), associated with a body A which is related to a global inertial frame 
through the transformation 

x' ^z'j,(t) + X'^, (47) 

z^(f) being the arbitrary motion of the local A-frame. The equations of motion with respect to the local A-frame reads (Damour, 
Soffel, Xu 1993): 

^ + ^.0, (48) 
dt dX'^ 

dipAv'J d ■ I dV\. 

, ^ + — (pAv' vi + fJ) = PA— , (49) 
dt dX{ '^'^ ^ ^ ' dX' 

A ^ 

where pA(f, Za) = p^it, Za) is the mass volumic density expressed in the local A-frame, being the velocity with respect to this 
frame while t'j denotes the stress tensor The following effective potential appears 

N j2 

V^^it, X^) = y y»(r, Za + Xa) - Vt,{t, za) - "t^ -^a , (50) 

^ df2 

where 

y4(f,ZA) = ^y^f,ZA), (51) 

B^A 

the considered body A being tidally interacting with A^- 1 perturbing extended bodies B; is the potential of each body B different 
from A. The last term of Eq. (50 1 re pres ents the inertial effects on the accelerated local frame A. This effective potential can be spht 
into the potential of A given in Eq. (35 i and a tidal potential, Vj, as follow: 

yA^yA + yA^ (52) 

the tidal part being given by 

y^(f, Xa) = y4(f, ZA + Xa) - Vt,{t, Za) - ^ ■ Za . (53) 



By integrating over the body A the equations of motion Eqs. (48 1 and (49 1, we get the equations for the conservation of the total 
mass of A, M'^ - Ma, and for the second time derivative of the local dipole moment M'j^ (Damour et al. 1992, Hartmann et al. 
1994): 

5^=0, ,54, 

at 



^^Pa^^/'Za. (55) 



dfi Ja dX-^ 



We can now derive the expression of Vj into series by using an STF expansion of ascending powers of Za as follow 

- E Z TTJ^aGI > (56) 

where G^" are the local effective tidal moments with 



Gf = d^V^JzA) - ^ , (57) 

Gti,-d!,..,,vUzA). (58) 

We assume that the origin of the local A-frame coincide with the center of mass of A, i.e. the dipole moment M^it) vanishes. Since 
now the {/, m] indices are related to their associated body: for example for A we use {Ia, '«a) while for B we use {Zb, '«b)- With these 
definitions, the right-hand side of Eq. (55 i can be written as: 

/a^^^''^^ = 2Z^<^1- (59) 

-'a SM/I^j'A- 

The local equation of motion is finally obtained from the d' Alembert Principle with 

^M,^0 = 0, (60) 
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which leads, by using Eqs. ( [59| and ( [57| , to: 



d/2 



/a! 



A 



(61) 



Expressing in term of Eq. ( 32 1 as in Hartmann, Soffel & Kioustelidis (1994), we finally obtain 

B^AbA=0/B=0 ^ \'AB/ 



(62) 



Using Eqs. ([57]i-([58]l and (62 1, the tidal potential can be expressed as: 



1 

/"AB 



(63) 



Using Eq. ( 15 1, we get 



1/a=2/b=0 '^AB 



,^ (2/a+2/b-1)! 



/a = 1 'b=0 



'a+'b + 1 
AB 



(64) 



that leads, taking into account Eq. ( 33 1, to: 



-f7VlVV V ( iyA (2^A+2ZB-l)!! B f/fB.>HBt^A "^A^Lb 



fi^A V/A=2/B=0mB^-/B 

oo oo 



-j^Z Z Z Z<-»" 



/a=1 mA=-lA /B=OmB=- 



(2/a+2/b- 1)! 
(2Za-1)!!(2Zb-1) 



— Y' V'A.mA {/'B.mB/j, 

„^^/A,mA^^/B,mB^A^LA '^Lb 



"LaLb 
'a+'b + 1 



(65) 



We now replace the Cartesian multipole moments by their spherical harmonics representation. First of all, using Eq. ( 26 1 and putting 

Za = r = {r, 9, tp), we get: 

I 



3 

'«A=-1 

4r'*7r/A! 



A - (2L+1)!! ^ ^/a,,«a(^'W. 



(66) 



(67) 



mA=-lA 



Then, the last term of Eq. (65 i involving the product of basis function j/ can be written by using Eq. (21 



as follow: 



"LaLb 
..'a+'r + I 



^'a.MIA^'b.'MB "i-Ai-B 



^A i-B Ja+'b + I 



' AB 



(2/a-1)!!(2/b-1)! 
(2/a+2/b-1)!! 



/A.niA {//A+'B.mA+niB a 

V/B,mB"^iAiB 



/a+'b + 1 



Noting that the STF derivative in the RHS of Eq. ([68]l can be split into two parts 



di 



'A+fe + I 

V'^AB / 



^'a+'b 



— 5,(«LaLb) + nLALB5; 



'a+'b + 1 

^'ab 



and using the relation given by Hartmann, Soffel & Kioustelidis (1994) 



. , /a + /b + 1 . 2/a + 2Zb + 1 . 



we get 



di 



rAB 



'a+'b + 1 



'"AB 



2/a + 2Zr + 1 



B 

'a+'b+2 

'ab 



-WiXaZb ■ 



(68) 



(69) 



(70) 



(71) 
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Inserting Eqs. (|68]l and ( fTT) into Eq. ( [65| , using Eq. ( [25| with Za- Zb = ''ab = (''ab, ^^ab, i/'ab), we finally get the expression of the 
tidal potential for \r\ < Ra 



l\ cc mu Alt, V t a \ 

V V V V r iV-^ Af« jA.'«Ay. rr? ,, /^+'B.'"A+'"B(t^AB.yAB) 

2-i 2-i 2-1 2-1^ ' 2L + 1 

/A=2mA=-/A 'B=0;B=-mB ''ab 



-]tJ^Z Z £ Z £(-1)'^"tK24-2Zb.i)m,^,,,M,1„,^^^^^^^^^^^^^^ 



1 oo 'a °o 'b 



mA=-l /;=i m'^=~t'^ ;B=OmB=-/B 



l.fflA _ 
^+/b,M1a+'«B 



/ . +/b + 1 ,in^ +mB +mA 



(^abj ^Pab) 



L+/B+2 



(72) 

The respective physical meanings of terms I and n are clearly identified. Term I corresponds to the gravitational interaction of 
B with A, while term II is the acceleration responsible for the movement of the center of mass of A. In the case of a ponctual mass 
perturber B, we recall that we get (see for example Melchior 1971): 



V^(t, r, Tab) = (t, r, r^s) - r, Tab) 



where 



Mb being its mass. 



\r-rAB\ 



and yl(f,r,rAB) = G — 

'"AB 



1 + 



rAB ■ r 



'AB / 



Eq. (72 1 exactly corresponds to the Eq. (3.25) given in Hartmann, Soffel & Kioustelidis (1994) that we recall here: 



v,^(t,r,rAB) - Z Z Z (-i)'^^r;:;::<,.B'-'^^u(^'^) 



5';A+/B,mA+mB (^ABjVAb) 



/a=2 mA=-'A 'b=0 mB=-. 

Ia oo 'l 



'a+'b + 1 

'ab 



Ma 



Z Z Z Z (-i)'^r;::::MA„, MB,„,x'a,. 



/a=1 irtA^-U /B=OmB=-/B 



J'/A+'B.mA+'KB (^AB.V'Ab) 



'a+'b + 1 

'ab 



It is then recast in its general spectral form, using that Vj is real and expanding it in the spherical harmonics for \r\ < R/^ 



oo /a 



Vj it, r, Tab) = ^ G ^ ^ [Ai;/A,mA - '"ab) + ^i;/A,mA if' rAB)] r'" F;a,«,a V) > 



B#A ;a=Omja=-/a 

where the coefficients Ai-j^j„^ and Ai;/^ are respectively given by: 



(-1) 



Ia 



An 



2L + 



* /a+/b,"IA+'«B 



(Sab, tf ab) 



/B=OmB=-/B 



'a+'b + 1 

'ab 



and 



(73) 
(74) 



(75) 



(76) 



(77) 



1 471 

'm, 



y y y y (-i)'l-'(24 + 2/B + i)rf'"'^fM;^ VrV"* , 

4=1 m^=-4 'B=OmB=-/B 
(^AB, V'Ab) 



/' +/B + l,m'.+mB+"'A 



/ +/B+2 
'ab 



(78) 



The more general form of the tidal potential being derived, we now express Ai;/^ and Ai /^ as a function of the Keplerian 
orbit elements of the perturber B. 

Here, we take into account the relative inclinations of the spin of each body with respect to the orbital plane. It is then necessary to 
define three reference frames, represented on Fig. (|3]l, all centered on the center of mass of the considered body A, Oa- 

- an inertial frame : {Oa, Xr, Yr, Zr), time independent, with Zr in the direction of the total angular momentum of the whole 
system Lxotai - ^-Orbital + ^BodyA + Yjk ^BodyBj which is a first integral (we are studying here the two bodies interaction between 
A and each potential perturber with k e 11, A^J). 
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- an orbital frame Ho ■ {Oa, Xq, Yq, Zq). We define here three Euler angles to link this frame to : {Oa, Xr, Yr, Zr): 

- /b, the inclination of the orbital frame with respect to (Oa, Xr, Yr); 

- ojb, the argument of the pericenter; 

- fig, the longitude of the ascending node. 

Let us finally define the last three quantities associated to the elliptic elements of body B's center of mass: ab, the semi major 
axis, e^, the eccentricity and Mb, the mean anomaly with Mb ~ tist, being the mean motion. 

- a spin equatorial frame 'R-e-j '■ {(9a,Xe, Ye,Ze). This frame is rotating with the angular velocity, Qa. This frame is linked to 
■^R : {Oa,Xr, Yr, Zr) by three Euler angles: 

- sa, the obliquity, i.e. the inclination of the equatorial plane with respect to the reference plane (Oa, Xr, Yr); 

- ©A, the mean sideral angle where ©a - dQA/df. This is the angle between the minimal axis of inertia and the straight line 
due to the intersection of the planes [Oa, Xe, Ye) and (Oa, Xr, Yr). 

- <pA, the general precession angle. 




Fig. 3. Inertial Reference, Orbital and Equatorial rotating frames C^^r, "Rq and "Re j) and associated Euler's angles of orientation. 

The Kaula's transform is then used to explicitly express all the generic multipole expansion in spherical harmonics in term of 
keplerian elements. Using the results derived by Kaula (1962), the following identity is obtained: 



5^/,;,, (6'ab,i^ab) 1 



'ab 

where the kij coefficients are given by: 



tTjTjTj '"•>4'« ^^a) Fij^p (/b) G,,p,, (cb) exp [/1'/,„y,p,^l . 



*B j=-l p=0 q 



21+1(1- \J\)l 



' \ 47T (l + \j\)l- 

d'j^ (ea) is the obliquity function which is defined as follow for j > nr. 



d\^(sA) ^ (-ly-"" 



(i + jy.(i-jv. 

(I + my. (I -my. 



My) 



l-J 



(cos Sa) , 



(79) 



(80) 



(81) 



the p^"'^^ (x) being the Jacobi polynomials (cf. Abramowitz & Stegun, 1972). The value of the function for indices /' which do not 
verify j > m are deduced from: 

(82) 
(83) 



or from their symmetry properties: 



d'j^„(7r + SA) = (-iy-'d'_j„(sA) 
4,„ (£a) = (-iy-"'d'_j,-,„ (sa) = dlj (sa) . 
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On the other hand, one should note that: d'j^^^ (0) - 6j,„. 
The inclination function, Fijp (/b), is defined in a similar way: 





m 




2 


2 


(cos if 


2 


1 


-2 (cos (sin 


2 





V6(cos f)" (sin l)^ 


1 


1 


(cos f ) -3(cos|) (sin 


1 





- V6 cos £ (cos f ) (sin |) 








1 - 6 (cos 1^ (sin jj 



Table 1. Values of the obliquity function d'. (e) in the case where 1-2 and j > m obtained from Eq. (81 1 (adapted from Yoder, 
1995). 



An (l + j)l 



where 



I l,m 



21+ 1 



4n 



moreover, the following symmetry property is verified: 

Fi,-j.P (h) = 



21+1(1- jV. 
^ [(I - m)] (I + m)]]^- 



2' [(l-m) /2]\[(l + m) /2]l 



(l-m)' 



(-1)'-^ 

(l + jV- 



(84) 



(85) 



(86) 



The usual value of these functions are given in Tab. Q. 

The eccentricity functions G/^,^ (sb), are polynomial functions having for argument (see Kaula, 1962 and Laskar, 2005 for their 
detailed properties). Their values for usual sets {I, p,q} are given in Tab. ([3]l. In the case of weakly eccentric orbits, the summation 
over a small number of values for q is sufficient (q e [[-2, 2J). More details can be found in the appendix of Yoder (1995). 
Finally, the phase argument is given by: 



'i'l.mj,p,g = (/ - 2p + g) Mb + ^l.mj,p,g (cOb,^^, ©A, <*a) 



where 



This can be also written as: 



0) 



l,m.j,p,q 



(I - 2p) wb + j (f^B - "^a) - ("©A + (' - '") ; 



^l,m,j,p,q - O" l,m,p,q ("B, ^^a) t + (ffl,m.j,p,q (^^B , ^^B ' "^a) , 

where we have defined the tidal frequency: 



m,p,q 



(I — 2p + q)nB — otQa 



and 



^l,m,j,p,q = (^ - 2/7) CJb + j {^*B ~ ^a) + {l-m) \ 



(87) 

(88) 

(89) 
(90) 
(91) 



The Kaula's transform allows us to express each function of Tab, e. of (^ab, 6'ab, Vab), as a function of the Keplerian relative 
orbital elements of B in the A-frame. Applying Eq. (79 1 to Aj;/^ ,,,^ and An;;^,™^ respectively given in Eq. (77 i and in Eq. (78 1, we 
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I j 





1 

2 

1 
1 1 

1 2 

2 
2 1 
2 2 



(I) 



-|sin2/+i 
I sin^ / 

3 



J sin /(I + cos/) 
-| sin /cos/ 
-| sin/(l - cos/) 
|(1 +cos/)^ 
I sin^ / 



(1 -cos/)^ 



Table 2. Values of the inclination function F/j^p (I) in the case where 1 = 2. Values for j <0 can be deduced from Eq. ( 86 1 (adapted 
from Lambeck, 1980). 



/ 


P 


9 




P 


1 




2 





-2 


2 


2 


2 





2 





-1 


2 


2 


1 




2 








2 


2 





1 - |e2 + . . . 


2 





1 


2 


2 


-1 




2 





2 


2 


2 


-2 




2 


1 


-2 


2 


1 


2 


F + 


2 


1 


-1 


2 


1 


1 










2 


1 








Table 3. Values of the eccentricity function Gi^p^g(e) in the case where / = 2 (adapted from Lambeck, 1980). Note that several 
combinations {l,p,q} have the same Gi^p^g value. 



get: 



4 'b 
Ay,,., = (-l)'^^^(l-<J/„o)(l-5u,i);^ X t2\K«r.\M-i^K,..] 



^ ' ' lB=QmB=-lB 



^ /a+/b /a+'b 

X 



^ TjTj '^'a+'b J^^>A+mB (^a) ^'/a+Zb-Zp (^b) G;a+;b,P,9 (^b) exp [-/*l';^+ZB,,„^+mBj;p,9] 

(92) 



„'a+/b + 1 

"b ;=-(/a+;b) ;j=0 ? 



S. Mathis, C. Le Poncin-Lafitte: Tidal dynamics of extended bodies in planetary systems and multiple stars 



13 



and 



ll;/A,mA - 



1 4;r 



^'^^'Z Z Z Z (-i)'^"K-2ZB.i)r;:2iMA„,,iexp 

/'.=! m'=-l\ 'B=OmB=-/B 



'a."'a 



/,+/b + 1 'a+'b + 1 



/^+/B,mA+''* 'a+'b+2 

"^B '■=-(4+'b + 1) ■'='^ " 

X exp [— '*P;^+/B + l,m^+mB+mA,'',i,«] ' 



(93) 



where as in Eq. (43 1 M? = Mf " + mJ'' and = Mf* + M?"* . 

Like in Zahn (1966-1977), the tidal potential can be splitted into two components. The first one, V^.^ (r, Tab), is stationary (i.e. the 
tidal frequency vanishes: cr = 0). It corresponds to the axisymmetric permanent deformation induced by B. In the case of a ponctual 

mass perturber and of a system where all the spins are aligned, Zahn (1966-1977) shown that = "^^3(1 ~ ^b) r^PiicosO). 
Then, the second component is the time dependent part of the perturbation, V^j (t, r, tab), for which cr 0. 



2.6. The two bodies interaction potential 

The mutual gravitational interaction potentiaj^of two bodies A and B is defined as: 



VA-B(f,rAB)= r y"(f,r,rAB)dMA. 

JMa 

Following Hartmann, Soffel, KiousteUdis (1994), its expansion on STF-tensors is given by 



(-1)'« 



/a=0/b=0 



1 

'"AB 



(94) 



(95) 



Using once again Eqs. ( 18 1 and ( 33 1, we get: 

Ia 



/A=OmA=-/A 'B=OmB=-/B I 



J'/A+/B,'nA+™B (^AB,<^Ab) 
'a+'b + 1 

'ab 



(96) 



Finally, using the Kaula transformation given in Eqs. (79 1, (80 1 and (87 1 as previously done for Vj, Va-b is expressed as a function 
of the obliquity, sa, and of the Keplerian orbital elements orB: as, sb and Ib'. 



\Ia ^/a,"»a 
Ib ,»ib 



^A-B = Gi: z z K™x,.B(-i)'^r;i 

/a=0 mA=-/A 'b=0 mB=-/B 
J Ca+'b) /a+'b 

X z z z 

'<'/A+/B,i''^i',mA+mB (sa) (h) Gi^+i^^w,b (^b) exp [i^, 



(97) 



i'=-(/a+/b) h'=0 b 



This interaction potential contains all multipole-multipole couplings. It is used in section 3.2 to compute the disturbing function 
aimed to study the dynamics of an extended body in gravitational interaction with A. 



Since all type of gravitational potentials have been examined, we now study the dynamics of a system of extended bodies. 



3. Equations of motion 

3.1. External gravitational potential of a tidally perturbed body 

The goal of this section is to derive the external gravitational potential of a tidally perturbed extended body A by an extended body 
B. This potential is the sum of the structural self-gravitational potential of A, V^(t,r), and of y.^(f,r, Tab), the tidally induced 
gravitational potential corresponding to the response of A to the perturbing potential Vj (f, r, Tab): 

Vel^t (t, r, Tab) = (t, r) + (f, r, Tab) (98) 

^ the denomination of Va-b as a potential is not very pertinent since it has the dimension of the product of a mass by a potential. However, we 
keep it to stay coherent with Hartmann, Soffel, Kioustelidis (1994). 
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with the following definition for Vs : 



(99) 



'a=0 lK=-mK 



where M^^^^^^ are the multipole moments of A in the case where it is not tidally perturbed by any other body, in other words, in 

the case it is isolated. By definition the external gravitational potential is harmonic; therefore y^j(f,r, Tab) verifies the Laplace 
equation: 

v2y^'^^^(f,r,rAB) = if |r|>7?A, (100) 



that directly leads to the same equation for V!^ (t, r, tab): 



v2y^(f,r,rAB) = if \r\>R, 



(101) 



Following Lambeck (1980), Neron de Surgy (1996), Neron de Surgy & Laskar (1997) and Correia & Laskar (2003a, 2003b), we use 
the classical Love numbers, k^^, which allow us to characterize the response of the body A to the tidal perturbation. The boundary 

conditions for (f, r, Tab) are: 



V^{t, W^O,rAB) = 

(t, \r\ = Ra, /-ab) = S/a Via (^ I*"! = -^a, ^ab) ' 



(102) 



where V/^ is the /J spherical harmonic of Vj (t, r, Tab)- We also recall that using Eq. (76 1 V.^ (f, r, tab) has been expanded as follow 
for |r| < Ra: 

Vt (f. r, Tab) ^ G J] [Aij^ma (t, J^ab) + Ai,i^,„,, (f, Tab)] r'^ Y,,,„.^ (0, <p) . (103) 
Using the well-known properties of the Laplace's equation, we search the solution for when \r\ > R^ of the form: 



/a 

Z 

;a=0 mA=-/A 



„/a + 1 



Inserting Eqs. ( 103 i and ( 104 1 into ( 102 1, the final solution of Vj is then derived with: 



Ta 

'a,"1a 



where mJ'^'^ and mJ'''^ are given by: 

/a./WA /a. 01a O J 



'a,«ia 



'a.'iia 



(104) 



(105) 



(106) 



The response of the body A, which is described by the Love numbers, is the adiabatic one. However, it is well known that 
an elastic as well as a fluid body reacts to the tidal perturbation with a damping and a time delay which are due to the internal 
friction and diffusivities (in other words to the viscosity, v, and the thermal diffusivity, K, in a non-magnetic body). That allows us 
to transform the mechanical energy into thermal one which leads us to the dynamical evolution of the studied system (cf. Fig.Q). 
Therefore, we introduce a complex impedance, Zt^ (v, K; *?/,), with its associated argument, (v, K; *Pi): 



Zt^ (y, K; T^) = \Zj, (v, K; ^P^) | exp [idj, (v, K; ^PJ] 



(107) 



which describes this damping. We thus substitute ^^|Zt^| exp [/5ta] to A;^ in the Eq. (106i. L corresponds to the indices of the 
considered tidal Fourier's moda^ The different modellings that can be adopted for Zj^ and 5ta are described in Alexander (1973), 
Zahn (1977) and CoiTeia & Laskar (2003a). 



Using Eqs. (92 1 and (93 i, the expression of m7*'' and mJ^'^ are obtained: 



/a.bia 



M, 



Ta;I 

'a, '"A 



(-1)'' 



4-71 



21a 



^—^ktyi'^' (1 - 5,„o) (1 - 5/A,l) X I^TA;;A,mA.L,(v, K; U.,,.„rM,p,,)\\Ml,„Jjt 



/a, Ma 
ma 



^ /A+fB + l ^'*+'«--^''^/'"A+mB (*^a) F l^+I^J,p (/b) Gl^+l^^p,q (Cb) exp [/ ((JTA;/A,mA,/.l(^' 'i'l^+l^,m,,+mB,j,p,q) ' *I'/a 



+lB,mA+mBj,P,q ' 



•K,..b)1' 



(108) 



Note that each tidal Fourier's mode have its own dissipation rate as it as been shown by Zahn (1966-1977). 
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Ma 3 



^^i+Xv+me 4+/B+2 ''4+'B + l,r'^,^,„' +„B+„A (^a) ^'4+/b + 1,.,. (^b) G f^^i^^^^,^, (efi) 



1 /«+'b + 1 



'A""''"A""B 'A 

"b 



xexp 



' ^'Jta;;a,»!a,Li *^'/A+/B + l,'nA+"'B+mA.'-,s,«) *^/a+'b + 1.' 



:mj^+mB+mA,r,s,u l'^,m'^ Olvl, 



■B V 



(109) 



where M? = M?" + M^" FlandM^ = A/f * + M> 



/b.wib 'b.»'b /b,i"i 



/a,"1a 'a.hia 'a.iia ' 



As for V^, can be splitted into two components. The first one V^j (r, tab) is stationary. It corresponds to the permanent 
component j for which the tidal frequency (cr) vanishes. The second component ''ab) is the time-dependent one that 

corresponds to for which cr ?i 0. 

Finally the external potential of A is thus written in its more compact and general form for |r| > Ra- 



oo /a 
Ia=0 mA=-'A 



r'A+1 



where 



'a.ma /a,"Ia 'a,"Ia 'a.^a 'a.'Ma 'a,»!a 



(110) 



(111) 



3.2. Disturbing function 

The goal of this section is to derive the disturbing function, Ka-c, due to a tidally perturbed body A, acting on a body C of which 
dynamics is studied and which can be different from the perturber body B (see Fig. Q). 




Dissipative adjustment 



Fig. 4. Classical tidal dynamical system. The extended body B is tidally disturbing the extended body A which adjusts itself with a 
phase lag due to its internal friction processes. The dynamics of a third body C (different from B or not) is then studied. Q.a, "Bj 
ric are respectively the spin frequency of A, and the respective mean motions of B and C. 

First, the disturbing function is related to the mutual gravitational interaction potential (cf Tisserand, 1889-1891; Correia 2001) 
through: 

??A-c(f,rAc) = -T^VA-C; (112) 



The tidal multipole moments of B due to A can be derived using the same methodology and substituting A to B for the perturber and vice- versa. 
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the sign being due to the potentials convention adopted here. 

Using the definition of Va-c given in Eq. ([96]l, we deduce the explicit spectral expansion of '??a-c in the spherical harmonics: 



/A=OmA=-/A /c=Omc=-/c V '"aC j 

are respectively the mass multipole moments of the body A and of the body C, while tac, 6'ac and <pAC are the 
spherical coordinates of the center of mass of the body C in the A-frame (cf. Fig Q). Then, using the Kaula's transformation given 



in Eqs. (79 1, ( 80 1 and ( 87 1, "Ra-c is expressed as a function of the obliquity, sa, and of the Keplerian orbital elements of C: ac, ec 
and Iq '. 

= Z Z Z {<,«u<.c(-i)s;- 

/a=0 mA=-/A 'c=0 mc=-lc 

J (/a+'c) /a+'c "I 

X z^+fc + i Z Z y! KlA+lc,vd[^mA+mc (^a) f fA+/c,i'.vt' fe) Gi^+i^,h,j, (cc) exp [i'FfA+/c,mA+mc,v,vv,fc] k (114) 
'^C i'=-(/a+/c) ^=0 A j 

Here, three types of gravitational interaction are treated in our formalism (see also Eq. ([TTT])). To describe them, one has first to 
consider the two causes of the multipolar behaviour of the gravitational potential of a body. The first is due to its internal structure 
and dynamics. In the case of a solid body, it is due to its proper asymmetry while in the case of a fluid mass, the internal dynamical 
processes such as rotation or magnetic field will break the ideal spherical hydrostatic symmetry of the body. The second is the 
deformation of the body due to its response to the tidal perturbation exerted by the perturber(s). In the case studied here, it is the 
response of the body A to the perturbation exerted by B computed in the previous section. Therefore, we split here the ^-indexed 



mass multipole moments of each body as in Eq. (43 i: 



't.M't h.mk lk,tnk 

where Mf' is the self-structural contribution of the body while M^'' is the tidal one. 



Mf„, =Mf' +mJ' , (115) 

Ih.tnh h,jnj- >,.m>, ^ ' 



The three type of gravitational interaction are thus identified. The first is the interaction between the structural mass multipole 
moments of each body, M^^'^ J ^ ^ with k + one should note that {/a - 0, - 0)-{/c = 0, mc = 0) is the classical interaction 
between Ma and Mc, Mc being the mass of C. The second corresponds to the mixed interaction between the structural and the 

S T / 

tidal mass multipole moments, M^ M^ ' . The third is the interaction between the tidal mass multipole moments of each body, 

ill ' 
M/ M/ . Therefore, the disturbing function could be splitted into three terms: 

"Ra-C = '^A-C;S-S it, rAc) + '^A-C;S-T (?, ^ac) + '^A-C;T-T (t, rAc) , (116) 

where '??a-C;S-s is the disturbing function associated to the structure-structure interaction, Ka-c j-s is associated to the tide-structure 
interaction and ^^a-c t-t, is associated to the tide-tide interaction. 



Inserting Eq. ( 115 i into Eq. ( 1 14 1, the respective Fourier expansions of 'Ra-C;S-s, ^a-C;S-t and ^a-C;T-t are derived: 



^ oo Ic 

= y y y y iMf^ Mf- (-D'Ay;-™- 

Mc ^ ^ ^ ^ ^ 'c.mc ^ ^ 'k,' 

/a=0 mA=-lA /c=0 mc=-/c 



X 

a. 



H+h + l Z '^'A+/c,i''^[Cm!+mc (*^a) F ;a+/c,v,vv (^c) Gl^+l^^^^b (^c) eXp [iV lA+lc,mA+mc,v,w,h] \, (117) 



^A-cT-s = -— y y y y kMf* mj^ +m}'^ m^ )(-i)'*r^"'* 



'a <» '( 
Ia=0 'nA=-'A 'c=0 mc=-/c 

X 



f ^ KlA+lc,vd[^,mA+mc (^a) F Ia+Ic,v,w Uc) Gl^^l^^„^b (^c) eXp [i^ lA+lc,mA+mc,v,w,b\ \ (118) 



Ja+/l 

v,w,b 

and 

oo /a °° 'c 

^--^-^ = -]i^Z Z Z Z {<:.A<%c(-i)"r;:- 

/a=0 mA=-/A /c=0 mc=-/c 



X "^XTtTI Z '^'A+'c,i '^lv)t!+mc (*^a) F Ia+Ic,v,w ih) Gi^+i^^„^b (ec) exp [i^ iA+ic,mA+mc,v,w,b] y (119) 

'^'^ v,w,b I 
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This classification of the three different types of interaction allows to explicitly generalize the classical case where the only 
considered extended body is the tidally perturbed one. A, while B and C are considered as ponctual masses. In this case, the 
interaction are restricted to the classical gravitational interaction between M?* , mJ^ , Mr, the mass of B, and Mr- 

6 'a,»1a' 'a,'«a' > L. 

By now, to lighten the equations, the tidal multipole moments of C are ignored. In a practical case, they have to be derived using 
eqs. ( 108|109 1 and taken into account. 

The disturbing function "Ra-c is thus reduced to the two first interactions: the respective structural moments of body A and of 
body C and the structural moments of body C with the tidal moments of body A. The Fourier expansion of the disturbing function 
is thus given by: 

(f,rAc), (120) 

lA_,mp^,lc,mc,y',w,b 



where 
with 

^S>-S\lp^,mp^,lc,mc,v,w,b 



it, rAc) 



-— Mf* Mf^ (-l)'*y/ 

'A. "'A 'C."'C ^ ' 'k 



Ik ^,'A,mA 
'c,mc 



Jk+lc 



ht+T '^'A+'ci^'^ICml+mc (*^a) F l^+lc,v,w (/c) Gi^+i^^,,.^h (sc) exp [iV l^+lc,mk+mc,v,w,b\ 



v,w,b 



(121) 



(122) 



and 



G 



Ij^+l^ + l ^/A+/c,v'4!mI+mc (^a) F l^+lc,v,w (h) Gl^+l^^^.b i^c) exp [i^ l^+lcjnA+mc,VM\b] • 



v,w,h 



(123) 



Using Eq. ( 105 1, the Tide-Structure interaction disturbing function is expanded as follow: 

'^T-S;/A,«iA,/c,mc,i',vv,fe - ^ (t, ''Ac) + ^ '^n;Li (t, ''ac) , 

lB,mB.j.p.q 



(124) 



where "Ri^ and 'Ri^l, respectively correspond to the mJ^^^^ and mJ^^^^ contribution to the external gravitational potential of the 
tidally perturbed body A, Vj. Inserting the explicit expansion for and Mj *^ given in Eqs. ( 108 i and ( 109 1 into ( 123 i, we get 



G 4-7T 



^2^771^''^^^'"' " ~ "^'-'^ |ZTA;/A.«AX,(v,^;^/A.fe,mA.,«Bj,p,.)lr;::::iM«,„j|M,^';,„^^ 



^Sc i„/A,niA 

:,mc 



'<^/a+/b j'^/A+'c.v'^;>t!+mB '^^a) d[!^mA-i-mc (^a) F (4) F l^+lc,v,w (^c) Glk+ls.P.g (^b) GlA+lc,n;b (^c) 



j'a+'c 



'a+'b + 1 'a+'c + 1 



Xexp [i[^U+lc,mA+mc,v,w,b ' *P/A+/B,mA+mBj>,9 + ^^fc,mf: ~ '^^/Ima + ^TA;lA.ir:A,L, {v, K^'V i^+i^^„^+,„gj^p^,^^j , 

where Li = {/b, ma, J, p, q], and for "Ri-Lt: 



(125) 



?Ad 



^n;Ln - ZTA;/A,mA,in(^'-^;*i'/^+/B + I,n4+mB+mA,r,i,»)-^il;Ln 



McMa 3 



I 

b.w'b' 



|mSc I /A,mA_ 
I Iciitc" Ic 



l'+lB,ni'+ma lc,mc" Icnic +/b+2 „'a+'c+1 'a 
flg "c 

X-^4+ZB+i,r,i (^b) F i^+ic,v,w (Ic) G,' . .5_„ (en) Gi^+i^^„m (ec) 



X exp 



' \ ^lA+lc,mA+mc,v,w,b ^ f +1^ + 1,1,1' +mB+mA 



.r.s.u+SM^^,,n^-6Mf^ ^,-6Ml „,^ +^Ta;/a,».a.l. (v, K' "i',;^,^^!,,,,;^™^^™,,,,,^) 



(126) 



where Lj = |/^, m^, Zb, "Jb, s, m}. ^^^J^^ and "^^^^^ correspond to the adiabatic response of A. Finally, as in Eq. J43jl M,^ 



M?" -hM^" andM^,,, = Mf* + M, 



Sa 



^Ta 



'b ,«ib 'b ,»ib 



/a,'«a 'a,"ia 'a,"ia' 
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3.3. Dynamical equations 

The external potential of the tidally perturbed body A by a body B being now well understood and known, our purpose here is to 
derive the dynamical equations for the evolution of the angular velocity (the angular momentum in term of Andoyer's variables) 
and the obliquity of A under the action of the gravitational interaction with a body C which could be different from the perturber 
B and of the Keplerian orbital elements of this third body: ac, ec and Ic. To achieve this, we follow the method adopted by Yoder 
(1995-1997) and Correia & Laskar (2003 a,b,c) who used the mutual interaction potential for the variation of the Andoyer's vari- 
ables and the disturbing function of the orbital elements. Here, gravitational interactions between B and C are not taken into account. 

Begining with the Andoyer's variables (cf. Andoyer, 1926), we respectively get the total angular momentum, L\ = Ia^a, Ia 
being the inertia momentum of A: 

dQA 

Ia^^Bq^Va-c, (127) 

and the obliquity, ba'- 



/a^^A^ cos fiA = -B^i^Va-C - COS SaBq^Va-C- 



(128) 



The classical equations of orbital evolution are given by the Lagrange's planetary equations (cf Brouwer & Clemence, 1961): 



at tiQUc 



(129) 



dec 
df 



rica^ec 



du)c^A-C + 



1 



-Bmc^A-C, 



die 



da-JRA-c + 



ncajec 

cos/c 



I ""c--"-^ I - 

^ 1 - sin Ic nca^Jl- sin Ic 



du>c^A-C- 



The Fourier expansion of the disturbing function "Ra-c is then introduced: 



(130) 
(131) 



l^,m/^,lc,mQ,\\w,b 



(132) 



where, as in Eq. \l\A\ , %f^,mUc,mc,v,w,b is given by: 
G 



Mc 



Mf" M? f-iv^y- 



I A jA,mA 

'c,mc 



1 



X ,^+/^+i ^/A+/c.'-'^!v«!+mc (*^a) F (Ic) Gl^+I^_»,^h (ec) exp [/*P/A+/c,;nA+"ic,i',H',i] ■ 



(133) 



Therefore, the following formal equations are obtained for the dynamical evolution of respectively the angular velocity of A and its 
obliquity: 

d^A X"' 
I a -7- = Mc > ^ Re {/ (ifiA + mc) 



IaO.A ^ cos EA = -Mc ^ Re [i [V + (niA + mc) cos Sa] %A,mA,lc,mc,v,w,b] 

lA,mA,lc,'nc,v,w,b 



(134) 
(135) 



where R^, (z) is the real part of a complex number or function z, while those for the keplerian elements are given by: 



J_ dOc = A J_ y {, (l^ +l^-2w + b) 'Rl^,mAJc,mc,v,w.b} 

ac at nc ai, , f—^ . 

C lA,mAjc,>nc,v,w,b 



(136) 



1 dec 



1 1 - e?. 1 



ec dt nc ej, aj, , , , 



Re< 



(Ia + lc- 2w) 



^cy 



+ b 



lA,inA,lc,mc,v,w,b 



(137) 



and 



— cos /c = — 
dt «c 



I 7 q2 

-i/l — e^-, C lA,mAjc,mc,v,w,b 



R,{i[v-(lA + lc-2w)cosIc]'Ri„, 



,mA,lc,mcA\w,b | 



(138) 
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If the tidal multipole moments of C are ignored, the following system is obtained using Eqs. ( |121| l and ( |124| l due to the dependance 
of Ki-Lj and "Ri-Ls on ©a and (pA- 

^A-^ = Mc ^A' ('"A + mc) -Rs- 



.Ad 
in 



(139) 



/Af^A ^ COS fiA = ^ Re [v + (Wa + mc) COS Sa] '^S-S: 



!A,tnA,lc,tnc,v,w,b 



-i ^ [{v - j) + (mc - me) cos sa] '^I;L, - ' ^ [(v - '") + (mc - (mA + nisjj cos sa] Him 



■ 2 [(^"^A +COSeA50^)ZTA;/A,,«AX,]^^^L, + [(<9<?a + COS Ea^Q^) ZTA;/A,».AXn] 

'b ,'"b j',/',? /' ,m' ,/b .ffiB ,r,s,u 



Ad 
;i.n 



(140) 



and 



1 dflc _ 2 1 
flc df tiQ 



^ Re]/(ZA + /C-2W + Z7) 



C /A,mA,/c>'"c>i%n',^ 



'^S-S;/A,mA,'c,mc,i',H'.i + ^ '^I;Li + ^ '^I;Z.n 



+ ^ ^Mc [ZTA;/A,mA.Li]'^I^L, + "^Mc [ZTA;/A,mAXi] "^n;! 

'b,'"bJ',/'.? 



(141) 



1 dec 1 1 - ec 1 



ec df «c ec '^c 



(Ia + lc- 2w) 



1 



- Z 



^TA;/A,mAXi 



lB,mB,j,p,q ;^,m^,/B,mB,r,.!,« 



L,m.,lB,mB,r,s,u 



1 



c ; 











2TA;/A,»!A,i.n 


<1 


> 


(142) 



d 111 

— cos Ir 

df 



I 2 E 

Vie I 2 ^ 

-Jl — e^, C lA.mAJc,mc,VM',b 



i [v - (Za + /c - 2w) cos /c] 



lB,mB,i,P,q l'j^,m'j^,lB,mB,r,s,u 



+ Z [("^"c - COS Icd^c) 

] ^I^L, + Z t*-'^"': ~ ^c5„c) ZTA;/A,mA.Ln] 

'b ,n'B ./P,? /' ,m' ,/b ./hb .r,.?," 



Ad 

M I 



(143) 



where the explicit expression for '??s-S;/A,«iA,'c,mc,i',w,fc' "^iXi ^^'^ '^nx, have been derived in Eqs. ( 122\ , ( 125 1 and ( 126 1 



4. Scaling laws in the case of an extended axisymmetric deformed perturber 

4.1. Comparison to the ponctual case 

Here, our goal is to quantify the term(s) of the disturbing function due to the non-ponctual behaviour of the perturber B and to 
compare it to the one in the ponctual mass case. 

To achieve this aim, some assumptions are assumed. First, we adopt the quadrupolar approximation for the response of A to the 
tidal excitation by B; thus, we put /a = 2 so that ^^nxn = 0. Then, we consider the simplified situation where the body of which 
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dynamics is studied is the tidal perturber: therefore B=C and Eq. ( |125| l becomes: 



G An 



B I 



X exp ; (5ta;2,hiaXi (v,/:;^'2+/B,mA+mBj,p,,)] • 



(144) 



On the other hand, since we are interested in the amplitude of we focus on its norm Finally, as we know that the 

dissipative part of the tide is very small compared to the adiabatic one (cf. Zahn 1966), we can assume that |Zt^| =s 1 in this first 
step. 

Let us first derive the term of \'Ri-^Li I due to the non-ponctual term of the gravific potential of B which has a non-zero average 
in time over an orbital period of B, (^n_p)j. ('') - ^/Tb ^ " p{t,r)dt that corresponds to the axisymmetric rotational and 
permanent tidal deformations (see Zahn 1977) (the same procedure can of course be applied to the non-stationnary and non- 
axisymmetric deformations, but we choose here to focus only on (^n-p)7- illustrate our purpose). Then, as the considered 
deformations of B are axisymmetric, we can expand them using the usual gravitational moments of B (Jig) as 



(r) = 



GMb 



+ (V^-P),^ where (v^-p),^ = G ^ {Mf^^ + mJ;,) ^^^f^ with Mf^^ + M,^^„ = - 



/b>0 



(145) 



Then, we obtain 

B flg 

(146) 

On the other hand, the term of associated to Mb, namely the disturbing function in the case where B is assumed to be a 

ponctual mass, is given by: 



b?^j (aB,eB,/B,eA): 



G 47r. 



M 5^X[r^;r] Kol -T[x2,]-[dl„,JeA)] [FupUb)] [G2,p,,(.b)] 
= ^ f ^4a^b ^ [xijf [dl„,, iSA)f [F2,J,p {h)f [G2,p,, (eB)f . 



(147) 



For this first evaluation of the ratio I'R''''^ I / |!??!'^'|, we focus on the configuration of minimum energy. In this state, the spins of A 



and B are aligned with the orbital one so that sa = /b = (that leads to j = ma and p - (2- ma + Ib) /2) and the orbit is circular 
(cb = 0). Then, we consider: 



I'Rll (flB,0,0,0)| 
teVflB,0,0,0)| 



(148) 



Using eqs. ( 146|147| i, we get its expression in function of 7;^ and of (T^b/ab) 

1 



4-7T 



1 ^2+;B,mA ^2+/b,2.* 



7V°y2.-A ^2,„,A i^2,2,0(0) 



Jf 

'B 



(149) 



As it has been emphasized by Zahn (1966-1977), the main mode of the dissipative tide ruling the secular evolution of the system is 
Ota - 2. We thus define £; such that 



/2.,«,2.!f(0) 



J. 



which can be recast into 



log 



3^2./.,2,^ (0) 



+ log Jib - Ib log 



Ob 
Rb 



Finally, keeping only into account the quadrupolar deformation of B (^2), we get: 



log (£2) = 2 



log I 2 I + log ■^2 - 2 log 



Ob 



(150) 



(151) 



(152) 



This gives us the order of magnitude of the terms due to the non-ponctual behaviour of B compared to the one obtained in the 
ponctual mass approximation. It is directly proportional to the squarred J2, thus increasing with (where en = Cl^/Q.^ with 
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= J^) in the case of the rotation-induced deformation and with Sj (where sj - q{R^la^)^ where q - Ma/Mb) in the 

tidal one, while it increases as (R^/a^)'^. Therefore, as it is shown in Fig. (jsjl, the non-ponctual terms have to be taken into account 
for strongly deformed perturbers (J2 > 10"^) in very close systems {ob/Rb < 5) while they decrease rapidly otherwise. This 
corresponds for example to the case of internal natural satellites of rapidly rotating giant planets as Jupiter and Saturn 
{J2 ~ 1.4697 10^ for Jupiter and J2 ~ 1.6332 ■ 10^ for Saturn; see Guillot 1999 and references therein). In the case of close 
Hot-Jupiters which are already synchronized (because of the tidal dissipation the rotation period is close to the orbital one), 
the rotation period is larger than 2 days (50 hours) that is roughly 5 times slower than Jupiter's rotation (10 hours). In this 
case, the flattening of Hot-Jupiter is less important and their J2 should be of the same order of the Earth's value (i.e. J2 runs 
from 10 to 10^). Then, the relative effect of the non-ponctual terms is less important. The situation may be different in the 
earUest evolutionary stages of those systems and this have to be studied in forthcoming works. 




2 4 6 8 10 



Fig. 5. Log £2 in function of ob/^b for J2 = 10'^ (red dashed Une), 10^ (purple long-dashed line), 10 ' (blue solid line). 
The non-ponctual terms have to be taken into account for strongly deformed perturbers (^2 > 10^) in very close systems 
(ob/^b ^ 5) while they decrease rapidly otherwise. 



4.2. Application to rapidly rotating binary stars 

Here, we apply the previous procedure to the case of rapidly rotating binary stars where en » fix- 
First, we roughly scale J2 of B to Jq as (cf Roxburgh 2001): 



■'2 



Mb 

Mo 



the deformation being thus mainly due to rotation. Then, Eq. ( 152 1 becomes 



log (£2) ^2 



log - +log/2" + 31og 



log 



Mb 

Mo 



•2 log 



Qo 



2 log 



Then, we introduce the stellar homology relations for main-sequence stars (cf. Kippenhahn & Weigert 1990): 



R 



Rr 



M 



Mr 



TT - TV where z\ = 



a + A-2 
a + 3A 



with A - 



dine 
dlnpjj 



and a - 



dine 
dlnT 



e being the nuclear energy production rate per unit mass inside the star. We finally obtain: 



log (£2)* 2 



log - + logy^" + (3zi - Dlog + 21og - 21og -^^ 



Mr 



Q.r 



Rf 



(153) 



(154) 



(155) 



(156) 



For Main-Sequence stars burning their hydrogen A ^ I while considering the pp chain leads to zi ~ 0.465 (M-b/Mq < 1.5) and the 
CNO cycle to zi ~ 0.795 (Mb /Mo > 1.5) (see Kippenhahn & Weigert 1990 for a detailed discussion). 

In Fig. (j6|, we plot £2 in function of Mq/Mq and Q.b/Q.q for ab - 3Rb taking 7® » 2.2 10"^. It is shown that the non- 
ponctual terms are not always a perturbation compared to the ponctual mass one (£2 > 10-^) for Mb > 40 Mo and Qb > 45 Qo 
that corresponds to rapidly rotating close binary massive stars. This has to be taken into account in future studies dedicated to the 
evolution of such stars. 

Since, we have now identified the regime where the non-ponctual terms have to be taken into account, those have to be examined 
by integrating the complete set of dynamical equations. This will be presented in a forthcoming paper. 
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Fig. 6. Log £2 for l^^o (red dashed line), 40Mq (purple long-dashed line) and 120Mq (solid blue line) Main-Sequence stars in 
function of ilfi/fi© at ob = 3/?b taking j!^ ^ 2.2 10^. Non-ponctual terms are not always a perturbation compared to the 
ponctual mass one (£2 > 10^) for Mb > 40 Mq and D.b > 45 that corresponds to rapidly rotating close binary massive 
stars. 

5. Conclusion 

This work represents another step towards the modelling of the tidal dynamical evolution of planetary systems and of multiple stars 
by taking into account the extended character of the bodies. We have used the STF tensors which allow us to treat in an analytical 
and compact way the complex couplings between the multipole behaviour of gravitational fields of two extended bodies. First, the 
main properties of STF tensors have been defined and derived. Next, after deriving the classical multipole expansion of the self 
gravitational field of an extended body A, we have derived the tidal potential related to the interaction between two extended bodies 
A and B, which is a generalization of the previous results where the perturber is treated through the ponctual mass approximation. 
Using those results, the external potential of a tidally-perturbed body has been provided and used to obtain the associated mutual 
interaction potential with a body C for which the dynamics is studied and which could be different from the previous tidal perturber 
B, as well as the corresponding disturbing function. The disturbing function being expanded into Fourier's series, the dynamical 
equations are derived. Those allows us to study the dynamical evolution of the angular momentum and the obliquity of the body 
A simultaneously with the Keplerian oribtal elements of C relative to the center of mass of A, namely the semi major axis, oq, 
the eccentricity, ec, and the orbital inclination, Iq- The equations have been derived in a general way without any linearization of 
function of any orbital elements as it has been done in Hut (1981-1982) and in Correia & Laskar (2003c). This allows us to study 
the strongly non-linear problem of inclined, eccentric orbits. Therefore, this formalism can be useful for elliptic systems. 

The major interest of this modelling concerns the dissipation of the tidal mechanical energy. Here, the response of each extended 
body is parametrized through the Love number ki for the adiabatic one, and an impedance Zj and its associated delay 6^ which 
describes the damping and the phase lag due to the (viscous and thermal) dissipation acting on the bulk induced by the perturber. 
Self-consistent modellings have to be developed for elasto-viscous bodies (Rogister and Rochester, 2004) as well as for fluid bodies 
(cf. Zahn, 1977) as well as equilibrium and stability conditions. This will be done for extended fluid bodies in a forthcoming sery 
of papers. 

Moreover, due to the general character of this work, various applications are possible. The first one deals with the dynamics of 
planets which are very close to their star as observed in several extrasolar planetary systems (cf. Laskar & Correia, 2004; Levrard et 
al., 2007; Fabrycky et al., 2007; Correia et al. 2008). Next, in our own Solar System, dynamical systems, such as giant planets and 
their internal satellites, can be studied by taking into account their spatial extension and their multipolar behaviour due for example 
to their rapid rotation. Finally, this work can be applied to the dynamical evolution of close binary stars where the tidal interaction 
and its dissipation dominate the behaviour of the system until it has reached its lower energy state, where the spins are all aligned, 
the orbits are circular, and the components are synchronized with the orbital motion. For these coming applications, the potential 
impact of the spatial extension of bodies on dynamics has to be carefully evaluated and understood. 
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